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Abstract 


A practical design tool which emphasizes the analysis of flowfields near 
the injector face of liquid rocket engines has been developed and used to 
simulate preliminary configurations of NASA's Fastrac and vortex engines. 
This computational design tool is sufficiently detailed to predict the 
interactive effects of injector element impingement angles and points and the 
momenta of the individual orifice flows and the combusting flow which results. 

In order to simulate a significant number of individual orifices, a 
homogeneous computational fluid dynamics model was developed. To describe 

sub- and supercritical liquid and vapor flows, the model utilized thermal and 
caloric equations of state which were valid over a wide range of pressures and 
temperatures. The model was constructed such that the local quality of the 
flow was determined directly. Since both the Fastrac and vortex engines 
utilize RP-1/L0X propellants, a simplified hydrocarbon combustion model was 
devised in order to accomplish three-dimensional, multiphase flow simulations. 

Such a model does not identify drops or their distribution, but it does allow 
the recirculating flow along the injector face and into the acoustic cavity 
and the film coolant flow to be accurately predicted. 


Introduction 

To provide detailed design information on liquid rocket engine injector 
performance, a computational simulation which relates the configuration of the 
injector faceplate to the flowfield generated in the engine is required. 
Injector elements are generally designed to form overlapping sprays to mix the 
fuel and oxidizer streams as they enter the combustor. Most engines operate 
at conditions such that the sprays are supercritical, or soon become so as 
they enter the combustor. Historically, patternator cold-flow experiments are 
used to estimate inlet mixture ratio distributions, ' and the JANNAF standard 
axisymmetric finite-rate (TDK) computer code is used to predict performance. 

Computational fluid dynamics (CFD) codes, like the FDNS-RFV code used 
herein, can be used to replace the TDK code without compromising the results. 

If the TDK code is used, the initial mixture ratio distributions must be 
estimated. The FDNS-RFV CFD code calculates these distributions directly. 
Even when using the CFD codes, three-dimensional effects are seldom simulated. 

Despite the major effort that has been devoted to analyzing spray effects by 
tracking drops, such analyses are too computationally intensive to describe 
the flows from the complex injector configurations in real engines. In 
addition, the underlying physics of stripping and atomizing supercritical 
liquid drops, where the surface tension does not exist, is not well 
understood. This work approximates the full spray dynamics by using a 
homogeneous spray model which retains the liquid-vapor flow features by using 
a real fluid thermodynamics model that describes the local quality of the 
fluid. Such an approximation is essential in order to simulate many of the 
important three-dimensional aspects of the combustor flow. The attributes o 
this model and examples of its use are presented herein. 
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Methodology 


To describe sub- and supercritical spray combustion with a homogeneous model, 
multiphase, multicomponent, reacting turbulent flows must be modeled. The methods of 
representing these effects are described as follows. 

Thermal and caloric equations of state, vapor pressure, heat of vaporization, surface 
tension, and transport properties are modeled with the equations of state proposed by 
Hirshfelder, et al 5 - 6 (we term these the HBMS equations of state) and with conventional 
correlations, 7 for the other properties. The property correlations used were not chosen for their 
absolute accuracy, but for their validity over a wide range of temperatures and pressures and for 
requiring a minimum of data to describe a particular species. These correlations are explicit in 
density and temperature. 


HBMS thermal equation of state: 
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HBMS caloric equation of state: 
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These equations are based on the "theorem of corresponding states" for real fluids, which 
essentially means that the p-v-T relations for all species are similar if these variables are 
normalized with their values at the critical point, i.e. if reduced values are used. The reduced 
values in these equations are indicated with a subscript r. H„ is the ideal gas species enthalpy. 
Z c is the compressibility for a given species at the critical point. The HBMS equations are 
attractive to use because arbitrary correlations for vapor pressure, heat of vaporization, and 
liquid densities can be used. Since multi -component fluid/ vapor mixtures may be present in the 
flowfield, the mixture properties are calculated by the additive volume method. This means that 
multiphase mixtures are treated as ideal solutions. To simulate RP-1, a surrogate fuel composed 
of several representative species was postulated. Properties of these pure species were then 
estimated and a mixture formed 8 . When not otherwise available, properties were estimated by 
using bond energies for their molecular structure with Benson’s methodology 9 . For conditions 
where the species become ideal gases, the thermodynamic data from the CEC code 10 or from the 
Thermodynamic Research Center at Texas A&M 11 were used. 


The CFD solver used was the Finite-Difference Navier-Stokes code with provision for 
using real fluid properties, the FDNS-RFV code. This code is pressure based; it differs from 
an ideal gas code in the methodology used to relate the pressure correction to the continuity 
equation and of course in the properties subroutines used. The pressure correction (p') equation 
used in the FDNS-RFV code is: 
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where the superscripts * and n denote the value at the intermediate and previous time steps, 
respectively. D p is the inverse of the matrix of the coefficients of the convective terms in the 
finite-difference form of the inviscid equations of motion. This is not an obvious definition, but 
is one which has made the FDNS-RFV code a useful solver. The sound speed used in the 
pressure correction equation is that calculated for the real fluid multi-component mixture. 

To describe both the Fastrac and vortex engines, a simplified RP-l/LOX combustion 
model was devised in order to accomplish three-dimensional, multiphase flow simulations. 
Global RP-1 pyrolysis, partial oxidation, and soot formation steps were used, and the 
combustion was completed with a wet CO elementary reaction mechanism. Hence, in addition 
to the hydrocarbons, the species CO, C0 2 , H 2 0, H, O, OH were considered. Using the wet CO 
mechanism allows for disassociation in the chamber and recombination in the nozzle, which is 
necessary to obtain accurate temperatures in the simulation. Averaged data from gas generator 
tests 12 and from subscale motor tests 13 were used to establish mean molecular weights, 
temperature, and soot levels for which the combustion model was tuned. The global submodel 
used for LOX/RP-1 combustion is shown in Table 1. The RP-1 is allowed to pyrolyze to C 2 H 4 
surrogate (not ethylene) which has the heat capacity of ethylene and the heat of formation of RP- 
1. Thus, the pyrolysis is thermally neutral, and the heat of combustion matches that of RP-1. 
The C 2 H 4 surrogate is used as the soot precursor. Soot is assumed to be C^H^, 
circumcircumcoronene, in order to treat it as a gas with a density like a solid. This model was 
designed to allow easier tuning of the reactions to match observed combustion product molecular 
weight and flame temperatures. The molecular weight and flame temperatures of a low mixture 
ratio used for validation were those observed by Aerojet in gas generator simulation 
experiments. Matching these properties also matches the sound speed of the combusting 
mixture. Soot oxidation was neglected in the simulations reported herein, because there was 
insufficient oxygen for the soot oxidation to occur since the global C 2 H 4 oxidation reaction was 
faster than the expected soot oxidation rate. A finite-rate vaporization model, such as that 
suggested by Johnson, et al 13 could be used to account for the heat and mass transfer effects 
between the RP-1 liquid drops and the surrounding hot gases, if sufficient data were available 
to justify its use. 

In all cases simulated, a k-e turbulence model was used to close the mass averaged 
transport equations solved by the code. Our experience is that this incompressible turbulence 
model overestimates the mixing in a combusting flowfield. However, since the liquid 
propellants are also mixed by this model, we concluded that there are currently insufficient data 
to better tune the turbulence model. Since a very fine grid is used to describe the details of the 
individual injector elements, the turbulence model essentially describes the subgrid turbulence 
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Table 1. Combustion Kinetic Models for RP-l/0 2 


RP-1 - 6.2 C 2 H 4 

Rate = AT B ex/>{-E/RT} [RP-1] 0 5 

A = 3.0117xl0 10 , B = 0, E/R = 2.523X10 4 


C 2 H 4 + 0 2 -* 2 CO + 2 H 2 

Rate = AT n exp{-E/RT} [C 2 H 4 ] 0 5 [OJ 10 

A = 1.29xl0 15 , B = 1, E/R = 2.516X10 4 


48 C 2 H 4 -* Soot + 84 H 2 

Rate = AT B exp{-E/RT} [C 2 H 4 ] 20 [0 2 + e] 05 

A = 5.1308xl0 12 , B = -2, E/R = 1.61 lxlO 4 , e = O.OOlp/Moj g-mole/cm 3 


effects. Such an analysis is a form of large scale eddy simulation, the large scale eddies being 
roughly of the size of the near flowfield of a single injector element. The homogeneous spray 
model was used to simulate a single element like-on-like (LOL) impinger injector element and 
a single element unlike impinger element for the configuration and flow conditions used in the 
cold-flow test in the JANNAF data base. Results for the LOL case are shown in Figure 1 , as 
well as the cold flow correlations. The mass flux distributions look very similar one-inch from 
the faceplate, and further downstream the correlations predict faster decay. The comparisons 
of mixture ratio distribution are shown in Figure 2. Since all of the injected flow was not 
recovered in the cold flow experiments and since the non-uniformity of the CFD predicted 
mixture ratio distributions appear more realistic, the CFD simulations appear to give as useful 
a result as the cold-flow data correlations. In all cases, regions of low mass flux should be 
considered immaterial. Furthermore, the cold-flow data are for a very limited range of 
configurations and flowrates. The CFD simulations could be used to parametrically and 
systematically evaluate the effects of impingement points, impingement angles, and injection 
velocities. Unless significant improvements are made in conducting the cold-flow experiments, 
such detailed information cannot be obtained. 

It could be argued that the homogeneous model does not accurately account for the 
mixing due to the interaction of drop clouds and that the liquid and vapor streams are not in 
thermal equilibrium. However, efforts to measure drop distributions in cold-flow LOL sprays 14 
have been unsuccessful because of the very dense drop concentrations in the spray proper. No 
spray combustion data exist with which to evaluate a model for the severe conditions which exist 
in a rocket motor. For typical rocket engines, the chamber conditions are mostly supercritical; 
thus, the surface tension of the liquid is reduced to zero and the heat transfer and drag 
coefficients between the drops and the surrounding gases becomes ill defined. For example, for 
LOX/RP-1 rocket engines, the chamber pressure is typically well above the critical pressure of 
RP-1, and for LOX the temperature quickly exceeds that of oxygen (frequently the pressure is 
also supercritical). Under these conditions the drops which do exist are far more unstable than 
those which are seen in cold-flow experiments. 15 Experiments conducted by the Air Force 
Research Laboratory, 16 show that supercritical injection of a round jet of liquid nitrogen into 
gaseous nitrogen behaves like a gas jet. This is precisely the behavior predicted by the 
homogeneous spray model. If viable data to describe droplet phenomena in high-pressure 
combusting flows do become available, the homogeneous spray model has two parameters (the 
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finite vaporization rate and the inter-element mixing rate) which may be adjusted to match such 
data. These changes would add negligible complexity to the CFD simulation. 


Numerical Simulations 

Two examples of the use of the homogeneous spray code will be presented: the Fastrac 
and the vortex motor. These engines were chosen because at the time this work was done they 
were under development by NASA/MSFC. The simulations were made for preliminary engine 
configurations to give insight as to how the injectors were expected to operate. The analyses 
presented do not represent the best and final design on either of these engines. 

The Fastrac Engine 

A segment of the Fastrac main combustion chamber near the comer of the faceplate and 
the acoustic cavity was analyzed. This segment simulates several like-on-like injector elements 
of the Rev.-C configuration and was chosen to evaluate the temperatures in the acoustic cavity 
and along the chamber wall. The number of elements were limited so that the simulation could 
be efficiently conducted on a PC. One injector element on row #8 and one element on row #7 
were simulated. The number of film coolant holes was arbitrarily increased by 6% to make the 
computational domain symmetric in the circumferential direction. The first 2.5 inches 
downstream of the faceplate were simulated. About 110K grid points were used for this 
simulation. This amounts to about 25 grid points per orifice opening. The boundary conditions 
and injector geometry and flow rates (per hole) are presented in Table 2. 


Table 2. The Inlet Flow Conditions and Configurations of the Fastrac Injector Rev.-C) 



Row tn 

Row #8 

Coolant 

Hole 

Oxidizer 

Fuel 

Oxidizer 

Fuel 

Pressure (psi) 

633 

633 

633 

633 

633 

Temperature (°R) 

190 

540 

190 

540 

540 

Density (lbm/ft 3 

67.14 

46.44 

67.14 

46.44 

46.44 

Diameter (inch) 

0.105 

0.069 

0.105 

0.069 

0.042 

Mass flow rate (lbm/s) 

0.3839 

0.1381 

0.3965 

0.1381 

0.0625 

Impinging angle 

O 

o 

O 

O 

40° 

o 

O 

uo 

45° swirl 

Cant angle 

8° 

-25° 

-16.5° 

-25° 

18° 


The predicted velocity and temperature fields are shown in Figure 3. The mean 
temperature in the cavity is 2000°R which is in good agreement with measured values. The 
isotherms are essentially the mixture ratios in this region. Notice the numerous small eddies and 
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recirculation zones which predicted by the analysis; these constitute the 
"large scale" eddies which are a result of the numerical solution of the 
transport equations. As the liquid propellants jet into the chamber the 
streams collide and mix. The combustion occurs where the propellants have 
intermingled enough to form a combustible mixture. As the combustion occurs, 
the hot gases expand and offer additional resistance to forming a well mixed 
flow. The simulation indicates extreme unmixedness, which is probably the 
case near the faceplate. Note that the cross-winds near the injector face are 
predicted, and these velocities create the flow which is induced into the 
acoustic cavity. The solution procedure was initialized by a cold flow 
calculation performed by not allowing the flow to react. The flame initiation 
was accomplished by assigning a hot spot (a few grid points were set to 
1500°K) near the wall. The major advantage offered by this simulation is that 
the injected propellants enter at the correct velocity, momentum, and energy 
levels for the engine. (Therefore, the predicted temperatures, mixture 
ratios, and velocities should be realistic.) The downstream boundary 
conditions are determined by mass conservation. The side of the control volume 
adjacent to the flows from the elements closer to the centerline are treated 
with slip boundary conditions. The sides of the control volume bounding the 
circumferential direction are treated with periodic boundary conditions in 
order to account for the swirl introduced by the film cooling flow. Ideally, 
a pie - shaped section of the motor would be analyzed, but the configuration of 
the injector face was symmetric only in a 90° quadrant. Also, too many rows 
of elements were used for all of them to be simulated on a PC. Nevertheless, 
the computational domain used was sufficient to allow a valid simulation of 
the flow through the acoustic cavity and near the chamber wall. 

The Vortex Engine 

The concept of a vortex engine is to use the vorticity of the injected 
propellants to provide longer residence times in the chamber and, hopefully, 
better mixing. A preliminary design of such a vortex motor was analyzed. 
This motor introduced the propellants through 16 pairs of unlike injector 
elements, tangential to and around the circumference of the combustion 
chamber. The injected propellants therefore produced swirl in the chamber to 
enhance the residence time in the short chamber. The propellants are injected 
at a high enough velocity to create the swirl, and they do mix very well at 
the injection point. In a conventional motor, the many injector elements 
provide well mixed propellants, whereas the vortex motor depends on passive 
mixing of the propellant streams. A one-sixteenth pie-shaped section of the 
chamber was simulated. Periodic boundary conditions were used on the sides of 
the section. The simulation was conducted only part way down the expansion 
nozzle, but this was far enough to evaluate the mixing efficiency in the 
chamber. Operating conditions for this engine are expected to be a mixture 
ratio of 2.6, 4.15 lb m /sec total propellant flow, 1000 psia chamber pressure. 

The chamber diameter is 2.35 in, throat diameter 0.96 in, and exit diameter 
of 3.03 in for an expansion ratio of 9.96. The diameters of the RP-1 injector 
and the LOX injector are 0.04 and 0.06 in, respectively. The impingement 
angle is 36E between the unlike injector elements. The computational domain 
consists of 90K grid points. 

The predicted temperature field and species concentrations are shown in 
Figure 4. The mixture ratio contours would correspond exactly to the 

isotherms. Hence, the propellants remain largely unmixed for these flow 
conditions. Streamlines from the fuel side of the injector element (i.e. the 
orifices nearer to the faceplate), the oxidizer side of the injector element, 
and from a near-wall point located in between the fuel and oxidizer elements 
are shown in Figure 5. The 



simulation shows the flow to be very unmixed and to possess moderate swirl. The turbulence 
mixing was simulated with the incompressible k-e model, which is believed to overestimate the 
mixing. However, even with the somewhat longer residence time due to the swirl, efficient 
mixing is not accomplished with this design. The predicted thrust from the CFD solution is 80% 
of that obtained from a one-dimensional, fully mixed, equilibrium solution. The predicted mean 
chamber pressure was 940 psia, with an 80 psia pressure difference between the injector wall 
and the engine centerline. NASA is building and testing this preliminary design as a baseline 
for the vortex concept, so a final verification of the CFD predictions will be made. 

Closure 


A homogeneous spray combustion model for analyzing liquid rocket engines has been 
developed. The unique feature of the CFD model is that real configurations of liquid rocket 
engines can be simulated from the individual orifices in the injector all the way to the nozzle 
exit. Rational approximations are made to the spray dynamics involved to provide this 
engineering CFD model. Real fluid properties are used to represent the liquid, vapor, and two- 
phase flow regions in the engine. For these multi-component combusting flows, some of the 
components may be subcritical while others are supercritical. Examples of this methodology are 
given for a segment of the Fastrac combustion chamber and for a preliminary design of a vortex 
motor, both of which utilize LOX/RP-1 propellants. The examples shown were prepared on a 
single Pentium II, 450 mHz processor. Depending on the computer used a larger number of 
injector elements (or fraction of the engine) could be analyzed. However, even with the full- 
time use of a Cray a solution of an entire Fastrac engine is not yet possible. However, a 
complete analysis of the vortex engine with its 32 injector elements which are of good 
circumferential symmetry is now possible on a PC. The code reported herein is currently being 
modified to use parallel processing to expedite the analysis of more complex liquid rocket 
engines. 
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Top: Empirical Solution; Bottom: CFD Solution 


Figure 1 Mass Flux Distributions for Like-on-Like Impinger 
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Figure 2 O/F Ratio Distributions for Like-on-Like Impinger 








Fig. 3 Temperature Contours and Plane Tangent Velocity Vectors at Various Cross Planes 
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Fig. 4 Temperature, Velocity Vectors, and Species Concentrations in the Vortex Engine 




Figure 5 Selected Streamlines of the Vortex Engine 


